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ABSTRACT 



We measured the angular diameter of the exoplanet host star t Dra with Geor- 
gia State University's Center for High Angular Resolution Astronomy (CHARA) 
Array interferometer, and, using the star's parallax and photometry from the 
literature, calculated its physical radius and effective temp erature. We the n 



combi ned our results with stellar oscillation f reque ncies from IZechmeister et al. 



(120081 ) and orbital elements from iKane et al.l (120101 ) to determine the masses for 
the star and exoplanet. Our value for the central star's mass is 1.82±0.23 M , 
which means the exoplanet's minimum mass is 12.6±1.1 Mj up i ter . Using our new 
effective temperature, we recalculated the habitable zone for the system, though 
it is well outside the star-planet separation. 



Subject headings: infrared: stars, stars: fundamental parameters, techniques: 
interferometric, stars: individual: HD 137759 



1. Introduction 



Frink et al.l (120021 ) announced the discovery of a substellar companion to i Draconis 
(K2 III, HD 137759) with a period of 536 days and a minimum mass for the companion of 
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8.9 Mj U piter- They used a stellar mass of 1.05 M from lAUende Prieto fc Lambertl (119991 ). 
who compared the absolute visual magnitude an d (B — V) color from Hipparcos data with 
theoretical isochrones from iBertelli et al.l ( I1994J ). However, Frink et al. acknowledge that 
evolutionary tracks for a range of masses are close together on the H-R diagram, so any 
slight change in the evolutionary model can have a large impact on the derived mass. 



Zechmeister et al.l (120081 . hereafter Z08) observed l Dra in search of stellar oscillations 



using three separate instruments over almost 8 years in order to refine the orbital parameters 
of the planet and determine the mass of the central star. They found low amplitude, solar- 
like oscillations with a frequency of 3.8 d -1 in two of the datasets and derived a stellar mass 
of 2.2 M using the equations 







L/L 

M/Mq 



■ 0.234 m/s 



and 



fn 



62 d 



^(Teff/5777 K) 



3.5 



(2) 



u osc / m s 1 

where i/ osc is the oscillation velocity amplitude, / max is the frequency of the strongest mode, 
and T e ff is the effective temperature. They used a l umino sity of 64.2±2.1 L & from the Hip- 
parcos catalog and T e g = 4490 K from iMc William! (Il990l ). Z08 then compared their 2.2 M Q 
value to those derived using T eff , surface gra vities (log q), a nd m etalli cities ([Fe/H]) f rom t he 



literature and the PARAM stellar model by iGirardi et al.l (120001 ) and Ida Silva et al. 



The m asses ranged from 1.05±0.36 M bas ed on values from 



( 119991 ) to 1.71±0.38 M Q based on values from lSantos et al.l (120041 ) . Because Z08's mass was 



Allende Prieto fc Lambert 



significantly higher than those derived using the model, they chose a mass of 1.4 M & to 
calculate the minimum mass of the companion, which they list as 10.3 Mj up it r- 

A more accurate way to estimate the star's mass would be to investigate the frequency 
splitting (A/o) using the equation 



A/ 



M/M & 



■ 11.66 d" 



(3) 



combined with an interferometrically measured radius, but unfortunately Z08's dataset was 
not suitable for measuring Af . 

The advantage interferometry brings is the ability to directly measure the angular di- 
ameter of the star. Then the physical radius can be determined using the distance from 



1 Available online at |http://stev.oapd.inaf.it/cgi-bin/param_1.0 
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the parallax, and T e g can be calculated. We combine our results with those from stellar 
oscillation frequencies to more completely understand the system through a description of 
the central star's and exoplanet's masses and the extent of the habitable zone. Section 2 
details our observing procedure, Section 3 discusses how t Dra's angular diameter and T eff 
were determined, and Section 4 explores the physical implications of the new measurements. 



2. Interferometric observations 

Observations were obtained using the CHARA Array, a six element Y-shaped optica l- 



infrared interferometer located on Mount Wilson, California (Iten Brummelaar et al.ll2005l ). 
We used the CHARA Classic and CLIMB beam combiners in the K'-band (2.13 jttm) while 
visible wavelengths (470-800 nm) were used for tracking and tip/tilt correcti ons. The observ 



i ng pr ocedure and data reduction process employed here are described in iMcAlister et al. 



(120051 ) . We observed l Dra over four nights spanning four years with two baselines: in 2007 
and 2008, we used the longest telescope pair Sl-El with a maximum baseline of 331 m and 
in 2011, we used a shorter telescope pair W1-W2 with a maximum baseline of 108 mj^| 

Chosing proper calibrator stars is vital because they act as the standard against which 
the scientific target is measured. We selected four calibrators (HD 128998, HD 139778, HD 
141472, and HD 145454) because they are single stars with expected visibility amplitudes 
>80% so they were very nearly unresolved on the baseline used. This meant uncertainties 
in the calibrators' diameters did not affect the target's diameter calculation as much as if 
the calibrator stars had a significant angular size on the sky. We interleaved calibrator and 
target star observations so that every target was flanked by calibrator observations made 
as close in time as possible, which allowed us to convert instrumental target and calibrator 
visibilities to calibrated visibilities for the target. 

To check for possible unseen close companions that would contaminate our observations, 
we created spectral energy distribution (SED) fits based on published UBVRIJHK photo- 
metric values obtained from the literature for each calibrator to establish diameter estimates. 
We combined the photometry with Kurucz model atmosphere^ based on T e g and log g values 
to calculate angular diameters for the calibrators. The stellar models were fit to observed 



2 The three arms of the CHARA Array are denoted by their cardinal directions: "S" , "E" , and "W" are 
south, east, and west, respectively. Each arm bears two telescopes, numbered "1" for the telescope farthest 
from the beam combining laboratory and 
distance between the telescopes. 



'2" for the telescope closer to the lab. The "baseline" is the 



'Available to download at http://kurucz.cfa.harvard.edu 
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photometry a fter c onverting; magnitudes to fluxes using IColina et al.l (119961 UBVRI) and 
Cohen et al.l (120031 JHK). The photometry, T eff and log g values, and resulting angular 



diameters for the calibrators are listed in Table [TJ There were no hints of excess emission 
associated with a low-mass stellar companion or circumstellar disk in the calibrators' SED 
fits (see Figured]). 



3. Determination of angular diameter and T e g- 



The observed quantity of an interferometer is defined as the visibility (V), which is fit 
with a model of a uniformly-illuminated disk (UD) that represents the observed face of the 
star. Diameter fits to V were based upon the UD approximation given by V = 2Ji(x)/x, 
where J\ is the first-order Bessel function and x = tiBOubX^ 1 , where B is the projected 
baseline at the star's position, 6*ud is the appar ent UD angular diamet er of the star, and 
A is the effective wavelength of the observation (jShao fc Colavital 119921 ) . Table [2] lists the 
Modified Julian Date (MJD), projected baseline at the time of observation (B), projected 
baseline position angle (0), calibrated visibility (V"), and error in V (crV) for i Dra. 

Figure [2] shows the UD fit to the observed visibilities and it is clear that this is not a 
sufficient model. A more realistic model of a star's disk involves limb-darkening (LD), par- 
ticularly in this case because most of the observations are on the third lobe of the visibility 
curve where secondary effects such as limb darkening play a m ore important role than in 
the curve between B=0 and the first null. lLacour et al.l (120081 ) analyzed multiple LD pre- 
scriptions for the K1.5 III star Arcturus and found the power law model to be a sufficient 
approximation. We chose to use this m odel becau s e Arc turus' spectral type is very close to 
that of i Dra. The model was based on iHestrofferl (119971 ): 



(4) 



where /(/i) is the brightness of a point source at wavelength A, 1(1) is the brightness at the 
center, /z = ->/l — (2r/# L D) 2 , r is the angular distance from the star's center, and 6*ld is the 
limb-darkened angular diameter. In terms of the star's visibilities, the power law prescription 
becomes 



fc>0 



r(a/2 + £; + 2)r(£; + l) 



where y r is the radial spa tial frequency, T is the Euler function (T(/c+l) = k\), a = 0.258 ± 
0.003 ( lLacour et al.l 120081 ). and k is the number of terms used. We tried a range of k values 
to check its effect on #ld and found that after 10 terms, #ld remained steady. Therefore we 
used k — 11. 
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The resulting UD and LD angular diameters and their errors (< 1%) are listed in Table 
[3j Additionally, the combination of t he interferometric measurement of the star's angular 
diameter plus the Hipparcos parallax (Ivan Leeuwenl 120071 ) allowed us to determine the star's 
physical radius, which is also listed in Table El 



For the #t,d fit, the errors were derived via the reduced x 2 minimization method (IWall fc Jenkins 



20031 ; iPress et al.lll992l ): the diameter fit with the lowest x 2 was found and the corresponding 



diameter was the final #ld for the star. The errors were calculated by finding the diameter 
at x 2 + 1 on either side of the minimum x 2 an d determining the difference between the x 2 
diameter and x 2 + 1 diameter. Figure [3] shows the LD diameter fit for i Dra. 

Once 6>ld was determined interferometrically, the T eS was calculated using the relation 



1 



BOL 



(6) 



where F BOL is the bolo metric flux and a is the Stefan-Boltzmann constant. -Fbol was de- 
termined in the following w ay: i Dra's V and K magnitudes were dereddened using the 
extinctio n curve describ e d inlCardeili et al.l (119891 ) and its interstellar absorption (Ay) value 
was from iFamaey et al.l (120051 ) . The intrinsic broad-band color (V — K) was calculated and 
the bolometric correct ion (BC) was determ ined by interpolating between the [Fe/H] = 0.0 
and +0.2 tables from lAlonso et al.l ( 119991 ). They point out that in the range of 6000 K 
> Teg > 4000 K, their BC calibration is symmetrically distributed around a ±0.10 magni- 
tude band when compared to other calibrations, so we assigned the BC an error of 0.10. The 
bolometric flux (Fbol) was then determined by applying the BC and the T e g was calculated. 
The star's luminosity (L) was also calculated using the absolute V magnitude and the BC. 
See Table for a summary of these parameters. 



Results and discussion 



We estimated the limb-darkened angular diameter for i Dra using two additional meth- 



ods as a check for our measurement. First, we created an SED 



it for the star as described 



in Section 2, where UBV photometry was from lMermilliodl (11991 ) ; RI p hotometry was from 



Monet et al.l (120031 ) , and JHK photometry was from ICutri et al 



( 20 03|). Figure B] shows the 



resulting fit. Second, we used the relationship described in lBlackwell fc Lynas-Grayl (119941 ) 
between the (V—K) color, T e g, and 9. Our measured 6*ld is 3.596±0.015 mas, the SED fit es- 
timates 3.81±0.23 mas, and the color-temperature-diameter relationship produces 3.63±0.53 
mas. Because i Dra is so bright in the K-band (0.7 mag) and because 2MASS measurements 
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saturate at magnitudes brighter than ~3.5 even when usi ng the shortest exposure tirneR we 
used the K magnitude from the Two- Micron Sky Survey ( INeugebauer &: Leightonlll969l ) for 
the color-diameter determination. 

The main sources of errors for the three methods are uncertainties in visibilities for 
the interferometric measurement, uncertainties in the comparison between observed fluxes 
and the model fluxes for a given T e g- and log g for the SED estimate, and uncertainties in 
the parameters of the relation and the spread of stars around that relation for the color- 
temperature-diameter determination. All three diameters agree within their errors but our 
interferometric measurements provide an error approximately 15 and 35 times smaller than 
the other methods, respectively. 

With our newly calculated T eff , we were able to estimate the mass of the central star 
using Equations 1 and 2 and obtained a mass of 1.82 ± 0.23 M . We then cal c ulated the 
exoplanet's minimum mass using the orbital parameters presented in lKane et al.l (120101 ) and 
the equation 



/(m) 



m sin i 



P 



(KVl-e 



2\3 



(7) 



(M + m) 2 2nG 

where the period P was 510.72±0.07 days, the amplitude K was 306.0±3.8 m/s, and the 
eccentricity e was 0.713±0.008. Our calculation produced a minimum mass of 12.6±1.1 
Mjupiter, which converged in two iterations. 

Our Tpff of 4545±110 K is within the range listed in Z08 , whic h spans 4466±100 K 
( Ullende Prieto fc Lambert! Il999h to 4775±113 K jSantos et al1l2004J ). On the other hand, 
the stellar mass derived here is lower than that calculated in Z08 (2.2 M ) and slightly higher 
than the range presented in their paper ( 1.2-1.7 M^), thoug h it overlaps within errors with 
the mass derived from the T e s and log g in lSantos et al.l ( 120041 ) . Hopefully future observations 
of this star will determine the frequency splitting (see Equation 3), which will allow for the 
direct measurement of the star's mass when combined with the interferometrically measured 
radius. 



Using the following equations from I Jones et al.l (120061 ) , we were also able to calculate 
the size of the system's habitable zone: 



and 



S b ,i( T cs) = ( 4 - 190 x 10_8 T cs) ~ (2-139 x 10~ 4 T cS ) + 1.296 
S b>0 {T eS ) = (6.190 x 10~ 9 T c 2 ff ) - (1.319 x 10~ 5 T cff ) + 0.2341 



(9) 



4 Explanatory Supplement to the 2MASS All Sky Data Release and Extended Mission Products, 
http: / / www.ipac.caltech.edu/2mass /releases/allsky/doc/[ 
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where Sb,i(T c s) and S&, (T c ff) are the critical fluxes at the inner and outer boundaries in 
units of the solar constant. The inner and outer physical boundaries r i o in AU were then 
calculated using 



We obtained habitable zone boundaries of 6.8 AU and 13.5 AU. i Dra's planet has a 
semimajor axis of 1.34 AU (Z08), so there is no chance the planet orbits anywhere near the 
habitable zone. 

The CHARA Array is funded by the National Science Foundation through NSF grant 
AST-0606958 and by Georgia State University through the College of Arts and Sciences, 
and by the W.M. Keck Foundation. STR acknowledges partial support by NASA grant 
NNH09AK731. This research has made use of the SIMBAD database, operated at CDS, 
Strasbourg, France. This publication makes use of data products from the Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts and the Infrared 
Processing and Analysis Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science Foundation. 
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Table 1. Calibrator Information. 



Parameter 


HD 128998 


HD 139778 


HD 141472 


HD 145454 


U magnitude 


5.84 


N/A 


N/A 


5.35 


B magnitude 


5.84 


6.94 


7.31 


5.42 


V magnitude 


5.85 


5.87 


5.92 


5.44 


R magnitude 


5.82 


5.21 


5.09 


5.46 


I magnitude 


5.84 


4.65 


4.23 


5.50 


J magnitude 


5.76 


4.15 


3.55 


5.37 


H magnitude 


5.80 


3.56 


2.85 


5.43 


K magnitude 


5.76 


3.30 


2.49 


5.43 


Extinction Ay 


N/A 


0.13 


0.13 


N/A 


T eff (K) 


9395 


4525 


4200 


9772 


log g (cm s~ 2 ) 


4.14 


1.94 


1.94 


4.13 


6> UD (mas) 


0.23±0.01 


1.08±0.08 


1.48±0.08 


0.27±0.02 



Note. — The ph otome tric values are from the follo wing sources : UBV 



Mermilliodl (119911 ). RI 



Monet et all (120031) JHK - ICutri et all (120031 ). 



The Ay value s are f rom iFamaey et al.l (120051 ) . All T e ff and log g values 
are from ICoxl (120001 ) and are b ased on the star's spectral type, except 
for HD 145454, which is from lAUende Prieto fc Lambertl (119991 ) . The 
uniform-disk angular diameters (#ud) are the result of the SED fitting 
procedure described in Section 2. 
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Table 2. t Dra's Calibrated Visibilities. 



Calib b e 

HD MJD (m) (deg) V aV 



145454 



145454 



128998 



145454 



139778 



141472 



54358. 


188 


306. 


.28 


115 


,2 


0. 


.041 


0. 


003 


54358. 


198 


304. 


.57 


117. 


,5 





039 


0. 


.003 


54358. 


211 


301. 


,74 


120 


,7 


0. 


.033 





.002 


54643. 


239 


311 


,24 


253. 


,6 


0. 


038 


0. 


005 


54643. 


252 


312. 


.42 


256. 


,8 





.043 





.004 


54643. 


265 


313. 


.31 


256. 





0. 


.040 





.004 


54643, 


276 


313 


91 


262. 


,9 


0. 


047 


0. 


008 


54643. 


287 


314. 


.27 


265. 


,5 


0. 


053 





.010 


54643. 


298 


314. 


.47 


268. 


,2 


0. 


041 


0. 


.007 


54648. 


261 


313 


.81 


262 


3 


0. 


.044 





007 


54648. 


273 


314. 


.23 


265. 


,1 





.046 





.008 


54648. 


283 


314. 


.44 


267. 


,7 





.044 


0, 


.008 


54648. 


354 


311 


.66 


105 


3 


0. 


051 


0. 


009 


54648. 


366 


310 


44 


108 


,2 


0. 


042 


0. 


005 


54648. 


378 


308. 


.91 


111 


,1 





038 


0. 


.006 


54648. 


389 


307. 


.30 


113. 


,7 


0. 


.042 


0. 


005 


54648. 


400 


305 


43 


116 


,4 


0. 


038 


0. 


005 


54648. 


410 


303. 


,34 


118. 


9 


0. 


.035 





.005 


54648. 


421 


300. 


95 


121. 


,5 





.030 





.004 


55775. 


159 


106 


.94 


184 


,9 


0. 


372 


0. 


048 


55775. 


169 


106. 


39 


188. 


,2 





.377 


0. 


.045 


55775. 


180 


105. 


.70 


191 


,6 





.384 


0. 


.032 


55775. 


192 


104. 


.84 


195 


3 


0. 


.356 


0. 


051 


55775. 


169 


106. 


.39 


188. 


.2 


0. 


.381 


0. 


.058 


55775. 


180 


105. 


.70 


191. 


.6 


0. 


.360 


0. 


055 


55775. 


192 


104 


.84 


195 


3 


0. 


343 


0. 


056 



Note. — The projected baseline position angle (©) is 
measured eastward from north. 
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Table 3. l Dra Stellar Parameters. 



Parameter 


Value 


Reference 


[Fe/H] 


+0.14 


Averaged from Soubiran et al. (2010) 


V magnitude 


3.29 ± 0.02 


Mermilliod (1991) 


K magnitude 


0.72 ± 0.04 


Neueebauer & Leiehton (1969) 


A Y 


0.03 


Famaev et al. (2005) 


BC 


0.42 ± 0.10 


Alonso et al. (1999) 


Luminosity (L Q ) 


55.3 ± 5.3 


Calculated here 


^bol (10 -8 erg s _1 cm -2 ) 


183.8 ± 17.7 


Calculated here 


6» UD (mas) 


3.434 ± 0.012 (0.3%) 


Measured here 


6» LD (mas) 


3.596 ± 0.015 (0.4%) 


Measured here 


-^linear (-R©) 


11.99 ± 0.06 (0.5%) 


Calculated here 


T eff (K) 


4545 ± 110 (2%) 


Calculated here 



13 





Fig. 1. — SED fits for the calibrator stars. The diamonds are fluxes derived from 
UBVRIJHK photometry (left to right) and the solid lines are the Kurucz stellar mod- 
els of the stars with the best fit angular diameters. See Tabled] for the values used to create 
the fits. 
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-0.10 b ... i ... i ... i ... i ... i ... i ... i 
300 302 304 306 308 310 312 314 

BASELINE (meters) 



Fig. 2. — i Dra uniform disk diameter fit. The upper panel shows the full visibility curve 
and the bottom panel close-up of the second lobe. The solid line represents the theoretical 
visibility curve for a star with the best fit #ud> the dotted lines are the la error limits of 
the diameter fit, the filled circles are the calibrated visibilities, and the vertical lines are the 
measured errors. The UD model is clearly insufficent to fit the visibilities. 
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Fig. 3. — l Dra limb-darkened disk diameter fit. The upper panel shows the full visibility 
curve and the bottom panel close-up of the second lobe. The symbols are the same as in 
Figure EJ 
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Fig. 4. — l Dra SED fit. The diamonds are fluxes derived from UBVRIJHK photometry 
(left to right) and the solid line is the Kurucz stella r model of a star with T c g = 4466 K and 
log g = 2.24 from lAllende Prieto fc Lambertl (Il999l ). The errors for the UBV measurements 
were less than 1%, no errors were listed for RI, and the errors for JHK were 20 to 30%, 
which are not indicated on the plot. 



